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QJ Abstract 
CD 

<~| We propose three different lattice operators to measure the intrinsic width £j of the chromo- 

electric flux tube in pure lattice gauge theories. In order to test these proposals we evaluate them 
for SU(2) and Ising LGTs in (2+1) dimensions in the vicinity of the deconfmement transition. 
Using dimensional reduction, we could perform the calculation in the effective 2d spin model 
using standard S-matrix techniques. We consistently found the same result for the three lattice 
operators. This result can be expressed in terms of the finite temperature string tension as follows 
£j = 2o(T) an< ^ implies that the intrinsic width of the flux tube diverges as the deconfmement 
transition is approached. 



1 Introduction 



In these last years much interest has been attracted by the study of the flux tube thickness in the 
confining regime of Lattice Gauge Theories (LGTs). At low temperatures the square width of the 
flux tube is predicted to increase logarithmically with the interquark distance [lj. This prediction 
has been confirmed by numerical simulations in various pure gauge theories, first in abelian models 
[2H5] and more recently also in non-abelian LGTs [6j[7] (see also [8] for some early attempt). 

The situation changes drastically as the deconfinement temperature is approached from below. In 
fact it can be shown that in this regime the dependence of the square width on the interquark distance 
becomes linear with a proportionality constant which diverges as the deconfinement transition is 
approached [9, 10j. Also these predictions were nicely confirmed by numerical simulations [lOyilj. 

All the theoretical predictions mentioned above were obtained using as effective string theory the 
Nambu-Goto action. However we know that the Nambu Goto action is consistent at the quantum 
level only in 26 dimensions and as a consequence for any LGT in three or four dimensions the correct 
effective string action, whatever it is, should deviate from the Nambu-Goto one at some high enough 
order in the perturbative expansion. The consequences of these deviations can be seen, for instance, 
looking at the critical behaviour of the theory in the vicinity of the deconfinement transition. For 
those models which undergo a second order deconfinement transition the Nambu Goto action predicts 
mean field critical indices while the Svetitsky-Yaffe analysis [12] (confirmed by a host of numerical 
simulations) predicts non trivial universality classes for the different LGTs, depending on the center 
of the gauge group. 

These deviations from the Nambu-Goto action may be of two independent types: they can be due 
to irrelevant or marginal terms in the effective string action or they can be due to the coupling of the 
massless degrees of freedom of the effective string with massive (non-stringy) "intrinsic" excitations 
of the flux tube. Due to our incomplete understanding of how the effective string description emerges 
from QCD we have no precise description of these massive modes, of their dynamical origin and of 
their action, but we have a few hints which may help our intuition. The most important one is that 
these massive modes should manifest themselves as a sort of "intrinsic width" of the flux tube. The 
Nambu Goto action in fact describes a string of vanishing intrinsic width and the flux tube thickness 
discussed in PQ (which we shall call in the following "effective width" to avoid confusion) is entirely 
due to the quantum fluctuations. In this respect the intrinsic width, which we shall denote in the 
following as £/, can be viewed as the residual thickness of the flux tube when the interquark distance 
R is pushed down to the scale (typically R ~ l/^/a®) below which the effective string description 
does not hold any more and does not contribute to the flux tube thickness. 

We expect that the intrinsic width should not depend on the interquark distance R (in opposition 
to the effective width which instead increases with R) but should instead be a new fundamental scale 
of the theory. It is tempting to identify this intrinsic width with the size of the Nielsen-Olesen vortex 
line (or equivalently the London penetration length of the dual Abrikosov vortices^) in the framework 
of the dual superconductor models of confinement which however has been proved rigorously only for 
Abelian Higgs or Georgi-Glashow like models. We shall further comment on this issue in the following. 
Let us also mention among the attempts to characterize the intrinsic width the holographic model 
recently discussed in [II] . 

Thanks to the remarkable universality theorems proved in |15pi6] we know that the first few per- 
turbative orders (in the expansion in powers of 1/oqR 2 of effective string action) are universal and 
that corrections to the Nambu-Goto action may appear only at very high order in 1/ctqR 2 , thus it is 
well possible that the most important source of corrections is not due to higher order non-universal 
terms in the effective string action but to the coupling to the massive modes. Indeed, corrections 

x See [13] for un updated discussion of this issue. 



1 



to the excited string states in (3+1) dimensions were recently observed in [T7] and associated to a 
possible coupling to massive modeq 2 ]. Moreover it was recently observed that in the 3d Ising gauge 
model corrections to the effective string potential larger than those predicted by the universality 
theorems are present both in the torus (interface) [19] and in the cylinder (Polyakov loop correla- 
tors) |20j geometries. All these examples suggest that a more detailed study of these intrinsic massive 
corrections would be important non only from a theoretical point of view but also to improve our 
comparisons with numerical simulations. 

However, performing this analysis in the original (d+1) dimensional LGT turns out to be very 
demanding from a numerical point of view. Thus in this paper we decided to adopt an alternative 
route. Following the approach that we recently developed in |21tl22|. we studied the intrinsic width 
in the vicinity of the deconfinement transition (but still in the confining phase) looking at suitable 
correlators in a d dimensional spin model which, in the framework the well known Svetitsky-Yaffe 
analysis [12] , is known to represent an accurate effective description of the original LGT in the vicinity 
of the deconflnement transition. 



2 Three proposals for the evaluation of the intrinsic width in the 
finite temperature LGTs 

Finding a lattice observable to measure the the intrinsic width of the flux tube at finite temperature 
is a rather nontrivial task. 

At low temperature the simplest and most natural option would be to identify the intrinsic width 
with the constant term which appears in the function which describes the dependence of the flux tube 
thickness on the interquark distance R. Following [TJ we expect at low temperature the following 
behaviour for the effective width (in adimensional units) 

w 2 = ^zl log{R/Rc) (1) 

where we denote with eg the zero temperature string tension of the gauge model. This prediction 
was confirmed by high precision numerical simulations both in the 3d gauge Ising model [2] model 
and in the 3d SU(2) Yang Mills theory [6]. 

The scale R c , measured in units of 1/y^oo turned out to be almost the same in the two models: 
R c ^ = 0.337(18) in Ising and R c ^ = 0.364(3) for SU(2). 

Eq.fpTJ) can be rewritten in adimensional units as 

a w 2 = ^1 log(i^) + b (2) 

The fit of the numerical data to this law gives very good x 2 values down to values of the R ~ 1/ y / oo . 
Deviations below this threshold may be assumed to be due to the intrinsic string thickness that 
we are looking for. In this way we would obtain both for Ising and for SU(2) £ry / <7o ~ ~ 0.4 — 0.5. 
The major drawback of this qualitative estimate is that the value of b rather strongly depends on 
the assumptions on the shape of the flux tube. This is not a gaussian (see for instance [23,24J for 
a detailed study in the 3d Ising case) and there is not a common consensus on which should be 
the functional form. Indeed, as we shall see below the two issues of a non zero intrinsic width and 
non-gaussian shape of the flux tube are deeply interconnected. 

The situation becomes even worse if one is interested in a finite temperature estimate of the 
intrinsic width since in this regime the square width of the flux tube increases linearly instead of 

2 Notice however that no deviation was observed in a similar calculation in (2+1) dimensions |18j . 
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logarithmically [9] and the constant term gets contributions also from the Nambu-Goto effective 
string action thus making it impossible to disentangle the intrinsic width. 

It is thus important to find an alternative observable which could allow a direct estimate of £/. 

Three possible candidates are: 

a] The large distance transverse behaviour of the plaquette - Polyakov loop correlator expectation 

value. 

b] The second moment in the transverse direction of the completely connected four point function 

of two Polyakov loops and two plaquettes |25| . 

c] The large distance transverse behaviour of the completely connected four point function of two 

Polyakov loops and two plaquettes. 

In order to characterize these proposals let us first review the standard way in which the effective 
width of the flux tube is evaluated on the lattice. In a finite temperature setting the lattice operator 
which is used to evaluate the flux through a plaquette p of the lattice is: 

/ p P't TT\ 

(ef>(p;P,P')) = K {ppn) P/ ~(U P ) (3) 

where P, P' are two Polyakov loops separated by R lattice spacings and U p is the operator associated 
with the plaquette p. Let us choose p to be equidistant from the two Polyakov loops (i.e. in the space 
perpendicular to the plane defined by the two Polyakov loops and intersecting this plane exactly at 
R/2). Then the width of the flux tube w is defined as: 



w 2 (R,L) = ^ (4) 



where R denotes the distance between the two Polyakov loops Nt the extension of the lattice in the 
compactified time direction (i.e. the inverse temperature: T = 1/iVt) and h the displacement of p 
from the P P' plane. 

With this definition the first proposal mentioned above corresponds to studying the \h\ >> R 
limit of \<j>{h\ R, Nt) ). If in this limit the correlator shows an exponentially decaying behaviour then 



we may estimate £/ as ((j>(h;R,Nt)\ ~ exp(—\h\/£i). The rationale behind this proposal is that in 

this limit the effective width, (which in the vicinity of the deconfinement transition is proportional 
to vR) becomes negligible and the scale which drives the exponential decay must be the intrinsic 
width. This proposal is clearly inspired by the definition of the London penetration length (see for 
instance |13| and references therein) and requires, to be defined, that the flux tube should not have 
a gaussian shape (as instead predicted by the Nambu-Goto action). This is another way to say that 
we are looking to effects beyond the Nambu-Goto effective string description. 

The second proposal above corresponds to looking at the four point function (see figfl]) 

(Mp,p'; p, p')) = (pp' ] u p u' p ) (5) 



where ( ) c denotes the completely connected expectation value and, as above, the two Polyakov 
loops are located at a distance R while the two plaquettes are located at the opposite sides with 
respect to the P P' plane at a distance y from the plane (We refer to fig[Tjfor the geometrical setting. 
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Notice, to avoid confusion that in figQ]we set r = R/2). This quantity can be considered, in analogy 
to the flux <f> defined in eq.([3]), as a sort of "intrinsic flux" <j)j. In fact this quantity becomes a 5 
function centered in y = if we try to evaluate it in a stringy framework (and in particular assuming 
the Nambu-Goto action in the physical gauge) in which any expectation value is obtained as a sum 
over single valued surfaces which cannot simultaneously pass through P and P' unless P = P' . 

Thus the second moment of 4>i is likely to be an estimator of the intrinsic width (in exact analogy 
with the use of the second moment of </> as an estimator for the effective width w). 

Finally, the peculiar properties of the quantity cj)j suggest another interesting possibility to esti- 
mate £/. In fact, it is conceivable to expect that the intrinsic width should be related to the change 
in the value of the correlation length between two plaquettes due to the presence of a flux tube (mim- 
icked by the correlator of the two Polyakov loops between them). The simplest way to measures this 
quantity is to study the ratio 



Then, according to the above conjecture, for any value of R this ratio is expected to decay in the 
large y/R limit as ~ exp(— |y|/£r) where £j is the intrinsic width. 

The main problem is that these observables are rather difficult to estimate numerically. For 
this reason, following the approach that we developed in our previous papers [SUES], we decided to 
address first the problem in the framework of the Svetitsky-Yaffe conjecture [12] looking at suitable 
correlators in 2d spin models. 

More precisely we chose to study gauge models in (2+1) dimensions with a second order decon- 
finement transitions and with a gauge group whose center is Z%. The two simplest cases belonging 
to this class are the (2+1) dimensional gauge Ising model and the (2+1) SU(2) Yang-Mills theory. 
These models in the vicinity of the deconfinement transition belong to the same universality class of 
the 2d Ising model. Moreover the confining phase below the transition corresponds in this mapping to 
the energy perturbation of the 2d Ising universality class. In the following section we shall construct 
a dimensionally reduced projection of the above correlators. We shall first review the general formal- 
ism and recall how the standard effective width may be obtained in this framework. Then we shall 
address the three proposal mentioned above to evaluate £/. Taking advantage of the integrability of 
the energy perturbation of the 2d Ising field theory, we shall be able to evaluate exactly the large 
distance behaviour of the correlators involved in these observables, and to eventually extract the 
value of the intrinsic width. We shall show that, as expected, the intrinsic width does not depend on 
the interquark distance and that the three different ways that we propose to evaluate it consistently 
give the same answer. 

3 Dimensional reduction and the Svetitsky Yaffe conjecture. 

According to the Svetitsky-Yaffe conjecture [12], if we choose to study a (d+1) LGT (with gauge 
group G) with a second order deconfinement phase transition and if the gauge group G is such that 
the d dimensional spin model with (global) symmetry group the center of G also has a continuous 
symmetry breaking phase transition then the two critical points must belong to the same universality 
class and we can use the spin model as an effective theory description for the (d+1) dimensional LGT 
in the neighbourhood of the deconfinement transition. This is the case if we choose for instance the 
(2+1) SU(2) LGT or the (2+1) Ising gauge model. Both have a second order deconfinement phase 
transition and in both cases one can use the 2d spin Ising model as an effective theory in the 
neighbourhood of the deconfinement point. 




(6) 
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In this effective description the Polyakov loops of the LGT are mapped into the spins of the Ising 
model, the confining phase of the LGT into the high temperature phase of the spin model and the 
plaquette operator of the LGT is mapped into the energy operator of the spin Ising model. The 
reason why such an approach is particularly effective in the present setting, is the integrability of 
the Ising field theory in zero magnetic field. This property allows for the calculation of multipoint 
correlation functions in the regime of large distances via the spectral expansion over form factors. 
The latter can be calculated exactly, and used as building blocks to write down analytic expressions 
for the large distance behavior of various correlators. In the following we will be interested in the 
asymptotic behaviour of three- and four-point correlators involving spin and energy operators. We 
refer the reader to 1 21 , 22 , 26-281 for technical details on the calculations. 



3.1 The 3-point correlator (aea): the first proposal to estimate the intrinsic 
width. 

Following the above discussion the operator which measures the flux density in presence of the 
Polyakov loops pair in the LGT is mapped into the three points function (acre), (see [22J for a 
detailed discussion of this mapping). In the particular case of the 2d Ising model this correlator 
can be evaluated in the vicinity of the critical point using the spectral expansion, and (see next 
section for further details and [26] for a review) leading to the following expression for the "flux" 
distribution [22J 

S(R, y) = (are) ~ e "W^+* . (7) 

where y denotes the transverse direction, m is the mass of the 2d Ising model and a large mR limit 
is assumed. 

The effective width of the flux tube is then given by the ratio 

This ratio is easy to evaluate in the in the large mR limit (see [10022] (Notice, to avoid confusion, 
that in [22] we used the variable r = R/2 and that we evaluated the unnormalized width, i.e. only 
the numerator of eq.([S]))) 

leading to the following result: 

w 2 (R)~ 7- + .... (9) 

where the dots stay for terms constant or proportional to negative powers of R. 

The intrinsic width in which we are interested can be obtained looking at the large y limit of 
ec L-(0)- I n this limit we obtain 



(Fff2R 



S(R,y) ~ l 'Vr e-^. (10) 



from which we find £/ = l/2m. As mentioned in the previous section the reason for which we 
could find a non trivial result for is related to the fact that S(R, y) has not a purely gaussian 
behaviour in the y variable. 
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3.2 The 4-point correlator (eaea): the second proposal to estimate the intrinsic 
width 

We want to study the 4-point correlation function of the form 

{e(x 1 )a(x 2 )e(x 3 )a(xi)) (11) 

in the high temperature phase of the 2d Ising model in zero magnetic field. Since we are interested 
in its large distance behaviour, we may use also in this case the Form Factors approach. 

To evaluate the intrinsic width we are interested in particular in the connected correlation function 
for a generic rhombus of side L, as described in figure [TJ In such a particular case, it is enough to 
consider the first two contributions to the spectral expansion (we proceed in close analogy with the 
analysis of [21] and [22]). 



(eaea)? = / / ' ' 3 (-Ff ) F^{9 12 - i(3 + 2ivr) F|(0 23 + i(3) e - mLt - cosh ^ +cosh ^>~ 2m y coshy3 + 

/oo in jn 



-mL(cosh #i+...+cosh 64) 



(12) 



where we listed the leading and next-to-leading contributions to the spectral expansion. Inserting 
the explicit expressions for the form factors, we have 



(eaea)t = -^H- f°° Ml.. d9 3 sinh(0 12 /2 - %0) sinh(l/2(0 23 + i/3)) 

-mL(cosh 9 

(2vrmFf 



1 —00 

g— mL(cosh6i+cosh63)— 2mj/cosh62 _j_ 
?a\2 /"oo 

. , . Ml-- d6 A sinh(0 12 /2) cosh(l/2(0 34 - 2i/3)) • 
^' l^ 7r J j-00 

•tanh(0i 2 /2) tanh(l/2(0 13 + 2i0)) tanh(l/2(0 23 + 2i/3))e~ mL ( cosh «i+--Hwhfl 4 ) + ^ 
Using the following property in the second term 

tanh(6>i 2 /2) tanh(l/2(0i 3 + 2if3)) tanh(l/2(0 23 + 2i/3)) = 

= tanh(6»i 2 /2) - tanh(l/2(0i 3 + 2i/3)) + tanh(l/2(0 23 + 2ip)) (14) 
and after some long but straightforward calculations we can reduce the correlator to this form 

„-2mL 

(eaea)? = (mK 7 ) 2 \K (2my) - sm(3Ki(2my)} + 

rriL 

+- / M\pp(e)K (2mLcosh{8/2))+qp(e)K 1 (2mLcosh(e/2))} + ... 

2! mL J-oo 

(15) 



where 



„ sinh 2 12 sinh 2 9/2 - sin 2 (3 

pp{6) = cos/3- 



cosh 6/2 cosh 9/2 cos 0(1 + tanh 2 0/2 tan 2 /3) 

_ sin 2 g 1 

9/31 j cos/3 cosh 2 + sinh 2 0/2 tan 2 /3 1 ' 
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Figure 1: Rhombus configuration for the four-point function (eaea)^ , eq. (fT5 

and 

y 



L = yV 2 + y 2 ; /3 = arcsin =. (17) 

yV 2 + y 2 

Paralleling the usual treatment for the flux tube thickness, we can compute the quantity 



JZdyy 2 (eaea )0 



™j = poo , , ttt- (18) 



using the leading, large-r behaviour of (eaea)f. In the limit r/y — > oo, we have 



f° dy (ever)* = (mFf) 2 (-e" 2 "" + . . .) (20) 
,/_oo \m 4rar / 

2 



Substituting in the formula for wj, we obtain 



, _ JZdyyHeaea)? i ^ 1 



in full agreement with the result of the previous sections. 

3.3 Large-?/ behaviour: the third proposal to estimate the intrinsic width. 

Following the third proposal discussed above, we are interested in computing the leading behaviour 
of the ratio (ecrecr)f /(e(y)e(O)) c in the limit of large transverse coordinate y. Introducing the variable 
x = r/y, and taking the small x limit we have 



<e(y)e(0)> t 



1 + u... 



2 % /vr(my) 1 / 2 



e~ 2my + ... (22) 
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where 



(e(y)e(0)) c = m 2 [Kf(my) - K*(my)] = ^-L^e"*"" + . . . (23) 

from which we can read off that the exponential decay is ruled by £/ = (2m) -1 , which can be 
interpreted as the intrinsic width. 



4 Discussion. 



4.1 Temperature dependence of the intrinsic width 

We have seen in the previous section that the three approaches consistently give the same result 
= l/(2m). It is interesting to convert this result in the language of the original (2+1) dimensional 
LGT. 

Looking at the behaviour of the correlator of two Polyakov loops and comparing it with its spin 
model projection we immediately identify (see for instance |10| ) 



m = a(T)/T (24) 

and hence 

«' " 2^f> < 25 > 

where T = 1/Nt is the temperature and <r(T) is the temperature dependent string tension which 
vanishes at T = T c and has a critical behaviour given by 

<r(T) = a (l-?r) . (26) 

where v is (following the Svetitsky-Yaffe conjecture) the critical index of the 2d spin model (i.e. 
v = 1) and do is the usual string tension (which in this framework is the T — > limit of <r(T)). Thus 
our analysis shows that in the vicinity of the deconfmement transition the intrinsic width depends 
on the temperature and (as the effective width) diverges at the deconfinement transition. 

Eq (j25p is the main result of our paper and in the vicinity of the deconfinement point, it may be 
rewritten as: 

T ( T ' 



It is interesting to observe that, using the well known relationship between T c and <To proposed by 
Olesen in the framework of the Nambu-Goto effective string model [29] . which in (2+1) dimensions 
reads T c 2 = 3o"o/vr we obtain in the vicinity of the deconfinement point 

«ry^~^(i-!)"' (28) 

which looks like a finite T extension of the result £/(T = 0)y / oo ~ 0.4 — 0.5 quoted above for the 
zero temperature limit of the intrinsic width 
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4.2 Non-analytic terms in the spin-spin correlator 

We mentioned in the introduction that the massive intrinsic excitations of the flux tube should 
manifest themselves not only in the intrinsic width but also as massive corrections to the interquark 
potential. The approach that we discussed in the previous sections can also be used to address this 
issue by looking at the subleading exponential corrections in the spin-spin correlator which is the 2d 
analogous of the Polyakov loops correlator. Performing a large distance expansion of the spin-spin 
correlator we find: 



(a(r)a(O)) 



cr\2 g-mr 



2ir (mr) 1 / 2 



1 + 



1 e 



-2mr 



64tt (mr) 4 



(29) 



from which we can obtain the non-analytic contribution to the interquark potential (T is the tem- 
perature in the gauge theory) 



V(r,T) 



rp i e -2mr 

log(cr(r)cr(0)) =mr + — log mr - ——- -j + .... 

2 o47r (mr) 4 



(30) 



Looking at the exponential term in the correlator we may easily extract the mass of this intrinsic 
excitation which turns out to nicely agree with the previous estimates, £/ = (2m) -1 . It is interesting 
to notice that a similar result was recently discussed in [33] in the framework of a holographic model 
for the intrinsic width. 
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